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Abstract 

We study the dynamical stability of the BTZ black string against fermonic and gravitational 
perturbations. The BTZ black string is not always stable against these perturbations. There 
exist threshold values for m 2 related to the compactification of the extra dimension for fermonic 
perturbation, scalar part of the gravitational perturbation and the tensor perturbation, respectively. 
Above the threshold values, perturbations are stable; while below these thresholds, perturbations 
can be unstable. We find that this non-trivial stability behavior qualitatively agrees with that 



oo 
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predicted by a thermodynamical argument, showing that the BTZ black string phase is not the 



privileged stable phase. 
PACS numbers: 04.30.Nk, 04.70.Bw 



I. INTRODUCTION 

Stability of gravitational configurations has been an intriguing subject of discussions since long. More recently, the 
discussion has been focused on the higher-dimensional black brane solutions provided by string theories. Gregory and 
Laflamme have shown in their pioneering work that neutral black strings in more than four dimensions suffer from 
a long- wavelength instability (GL instability) [H, [|[ . In the last few years there have been a lot of progress achieved 
in understanding various aspects of the Gregory-Laflamme instability such as the GL instability with the presence of 
charge 0, and the GL instability in the dS/AdS spacetimesj!, H[. A review on the black string instability and whole 
ly-^ ' list of relevant references are presented in [fj . 
■^j- . Besides dynamical stability, there are also thermodynamical stability problems of black holes. However for usual 
black holes, such as Schwarzschild holes etc, they are thermodynamically unstable but dynamically stable against 
, perturbations of metric. The profound connection between dynamical stability and thermodynamical stability was 
observed in black branes in supergravity. Gubser and Mitra proposed a correlated stability conjecture (GM conjec- 
ture) arguing that gravitational systems with translational symmetry lack stability against small perturbations often 
also lack thermodynamical stability 0]. Supports on this conjecture have been presented in [1,0] j while a set of coun- 
terexamples has also been raised j8|. It is of interest to examine this link further in various configurations. Recently a 
non-trivial relation between dynamical and thermodynamical properties in the background of charged Kaluza-Klein 
black hole with squashed horizons was discussed in [9(. 

In this work we are going to study the classical stability of the BTZ black string solution obtained by starting with 
a metric describing an accelerating black hole in AdS^, 



ds 2 = 



A 2 (x - y) 



(i) 



P(y) = -A + ky 2 - 2fhAy 3 , Q(x) = 1 + kx 2 ~ 2fhAx 3 , (2) 
which satisfies the Einstein equation with a negative cosmological constant 

Rab = -(3/l 2 4 )g A B, (3) 



and Li = l/(A\/\ + 1), A being the acceleration parameter. The spacetime is sliced onto the branes of x/y — const., 
and the visible brane is situated at x = 0. There are two configurations of interest: 
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(i) The black string phase. To orient with m = 0, the BTZ black string in AdS 4 can be obtained through coordinate 
transformation [1 01 ] . which is described by 

ds 2 = l\dz 2 + l 2 a 2 (z) [-H 2 dt 2 + H- 2 dr 2 + r 2 dcp 2 ] , (4) 

where H 2 (r) — \r 2 — k and a(z) — -\/Acoshz being the warp factor. We have the interest on A > where locally the 
constant z slices have the geometry of AdS^. If k = +1, surfaces of constant z have the geometry of BTZ black hole 
fl"l1 | described by 

ds 2 = -(Xr 2 - k)dt 2 + (Ar 2 - k^dr 2 + r 2 d(p 2 , (5) 

which is the solution of i? M „ = —2\g flt/ . 

Here, Greek letters take the values 0, 1, 2 describing the brane and the capital Latin letters take the values 0, 1, 2, 3 
describing the bulk, where 3 indicates the extra dimension (i.e. a; = t, x 1 = r, x 2 = ip, x 3 — z). 

(ii) The black hole phase. Turn to the case with ra > 0, Q(x) has only one positive root x-i- Taking Atp = 
4ir/\Q'(x2)\, the conical singularity can be avoided. The black hole localized on a brane is described by [l0( 



ds2 = 4 
A 2 



-l Xr 2 -k-—)dt 2 +(\r 2 -k-—) 1 dr 2 +r 2 d^ 2 



(6) 



which is the BTZ-like black hole, with extra terms of the form 2faA/r coming from the four-dimensional nature of 
the black hole. If 2mA <C A -1 / 2 , the extra term is negligible outside the horizon, and so the exterior geometry is 
essentially identical to that of a BTZ black hole [Io| . 

A global thermodynamical argument using the entropy was presented to argue about the stability of these two 
configurations [Tol j . For the BTZ black string, it was shown that its thermodynamic mass M4 precisely agrees with 
that of the projected BTZ black hole M3. The four-dimensional entropy of the black string was also found in agreement 
with the three dimensional entropy 



q _ tt /2M 3 

where M3 is the projected mass of the BTZ black string on the visible brane. The black hole entropy was obtained 
in the form 

Sbh= G 4 A 2 X + 3P + 2P- (8) 

The auxiliary variable z is introduced by the identity 2mA = ^^_J^y,]£ z an d noting that z is a monotonic function 

of fhA, which leads to fa A — > 0, z k, 2fhA\f\\ fa A — > 00, z « A 1 / 3 . 

Equating the black hole and black string entropies, Eq.([7]) and Eq.©, and identifying the mass of the two phases by 

substituting M 3 — \\^z^+2zzy ' we a f° ur th order polynomial in terms of the auxiliary variable z{\— z 2 — z 3 ) = 
0. The root z — simply indicates that the mass and entropy of the localized BTZ black hole (when faA — > 0) 
match that of the black string, which shows that thermodynamically there is no privileged stable phase in these two 
configurations. When faA increases so that the auxiliary variable z reaches z x , which satisfies z 2 + z 3 , — A, the black 
hole and the black string entropies equal to each other again. Except these two points, the entropy comparison of 
black hole and black string phase is not illuminating in the auxiliary variable z. However expressing the entropies of 
the black hole and the black string phase in terms of their masses, it was shown in Fig.l of [l(| that, for small faA, 
the black string phase entropy exceeds that of the localized black hole, while for large fhA, the result goes in the 
contrary; the changeover point is at the root z x . 

In the present work, we will study the stability of the BTZ black string against fermionic and gravitational per- 
turbations and examine the relation between the dynamical stability against different fields' perturbations and the 
thermodynamical stability observed in [Iol |. Since we have, in this background, exact solutions for both the black 
strings in 3+1 dimensions as well as BTZ black holes on the brane in 2+1-dimensions, we can examine the stability 
of the black string in detail. 



II. FERMION PERTURBATION 



In this section we take a first look at the stability of the BTZ black string against fermionic perturbation. We 
assume that the Fermion field exists in the bulk. We will derive the fermionic perturbation equation and examine the 
mode which signals the stability of the black string. 
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A. Perturbation Equation 



For the massless Fermion field, it is described by the Dirac equation[lz 



(9) 



Here 7" is the Dirac matrix satisfying {^ a ,j b } = 2r/ ab , rf h is the Lorentz metric with the signature (—1,1,1,1). e„ 
is the tetrad so chosen that e^e b 3 r] ab = g AB , whose indices in small letters are raised or lowered by r\ a h and indices in 
capital letters are raised or lowered by qab ■ r 'a is the affinity of covariant derivative, defined by Ta = § [7°, 7 b ]ef ebB-A, 
£bB;A = ebB,A — ^BA e bE- In the spacetime background of Eq.((4]), we can rewrite Eq.© as 



7° 

Separating the wave function by 
and taking 



(2Ar 2 - k) 



-0„¥ + 7 3 a d. 



2r{Xr 2 - k) 
*(t, r, <p, z) = R(r)Z(z)e~ luJt+l ^, 



3d z a 
2a 



* = 0. 



R(r) = [r 2 (Ar 2 -fc)]- 1 / 4 *^), 
3d z a 



Z(z) — exp 



2a 



-dz 



dz 
a 



we can obtain the radial wave equation in the form 

- iwi? _1 7 $ + H^dr® + i^r^ 1 ^ 2 + m~f 3 )$ = 0. 



(10) 

(11) 

(12) 
(13) 

(14) 



Here, \x takes the value of (integer times 2~k) / (period of ip) as required by the periodical condition in tp. The value 
of m is determined by the compactification of the extra-dimension. 
Choosing 
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1,2,3, 



and letting $ = (<f>i, 2 ) T , (j> a = ((j> a +,<f>a-) T (a = 1, 2), we have 

iu)H~ x ^>x + Ho 1 d r <\>\ + i(fir~ 1 cr 2 + mer 3 )0i = 0, 
-iuH^fa + H(T l d r 4>2 + iinr^a 2 + mer 3 )0 2 = 0. 

We focus on the second equation and the first one can be solved easily by changing u 
the Pauli matrices, Eq. (fTT)) gives 



(—iuiH 1 + im)4>2+ + (Hd r — jir 



= 0, 



Defining Y± = ( nn±m) 



±1/4 



(-iojH- 1 - im)(j)2- + {Hd r + A*r _1 )(/>2 + = 0. 

i>2+ ± 02-) j we can change these equations into 
(dr T S)Y± ± imY T = 0, 



where 



S{r) 



y/ r 2 m 2 + /i s 
mrH 
imr 
H\J r 2 m 2 + fi 2 



dr, 



fimH 



2{r 2 m 2 + /1 2 ) 



Decoupling Eqs. (fl"§|) , we have 



d 2 Y ± + (m 2 + S 2 ± dfS)Y ± = 0. 



(15) 

(16) 
(17) 

-u>. Explicitly writing out 
(18) 

(19) 

(20) 
(21) 

(22) 
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Eq. (fT9]) and Eq. ([22| allow us to look at the asymptotic behavior of Y± at the horizon r = rjj = y/k/X and at 
spatial infinity r — > oo. In the vicinity of the horizon r — > m + 0, Y± ~ (r — r#) p± , and from Ea.(|22p we require 
p± to satisfy 2pj_ — p± + lu 2 /2Xk ± iu)/2\Xk = 0. Since there exist only ingoing waves near the horizon, we have 
p + = 1/2 — iuj/2y/Xk, and p_ = — iu)/2y/~Xk and thus in the limit r — > r# + 0, 

y + ~ (V _ fH )3-?^s ) y_ _ ( r _ rH y^m_ (23) 

For a consistency check this result also satisfies the coupled first-order equation (|19|) . Similarly we let r — > 00 and 
find asymptotic solutions satisfying Eq. (fP9|) and Eq. ([22|) 

r -> 00. (24) 




B. Quasinormal Frequencies of Fermionic Perturbation 

We now investigate the stability problem by solving Eq. (119p . which is similar to the quasinormal modes (QNMs) 
study (a review on the QNMs study can be found in [13l.fl4l|). 

In AdS spacetime, the QNMs were defined in [TtJ by only ingoing modes near the horizon while at the spatial 
infinity the wave function here equivalent to r~ 1 Y±, vanishes. This boundary condition allows only a discrete set 
of complex lu to exist. If the imaginary part of the QNM frequency is negative, we will have the decaying mode. 
However, if the imaginary part of QNM frequency is positive, we will have a growing mode, indicating that the 
spacetime is unstable against the field perturbation. 

To calculate the QNM frequencies, we need to solve Eq. (fl9|) and select only solutions which satisfy \& — * as 
r — > 00. For simplicity, we will let fi — 0, which allows to solve Eq. pTI)) analytically. Working with the decoupled 
equation Eq. (|22p , we can rewrite them as 

H 2 d 2 Y- + Xrd r Y- + (m 2 - H 2 + iuH~ 2 Xr) F_ = 0. (25) 

This equation is symmetric under the transformation of r — > ar, m — > m/a and lu — > w/a. In the following we will 
set A = 1. Taking x = (r — l)/(r + 1), Y_ = x~ w f 2 (x — l)~ m y(x), Eq. ipS)) becomes a hypergeometric equation (see 
also H|) 



z(l - x)pt + [c - (1 + a + b)x] -f--aby = 0, (26) 
ax ax 



where a — —m, b = 1/2 — ito — m, and c = 1/2 — iu>. We choose its solution to be simply 2 -Fi(a, b, c; x), and 
eliminate the other one x 1 ~ c 2-^1 (a — c+1,6— e+1,2 — c;x), so that Y_ contains only ingoing modes near the 
horizon embodied by the factor x~ %M l 2 . Therefore F_ = x~ lu) (x — l)~ m 2-^1 (a, b, c; x). Considering 2-Fi(a, b, c; x) — 
(1 — x) c ~ a ~ b 2-Fi(c — a, c — b, c; x) [l5j], we can rewrite the solution as 

K_ = x~ iu) / 2 (x - 1)-H 2 F 1 (a,b,c;x), (27) 



where a = —\m\ and b = 1/2 — iui — \m\ satisfying Re(c — a — b) — 2\m\ > 0. We need to substitute Eq. (f27f into the 
first equation of Eqs. p^|) to get the explicit expression of Y + . Using 4- 2^1(0, b, c; x) = — 2 Fi{a + F b + 1, c + 1; x) 
we have 

Y+ = -iex^-^^ix-l)-^ \ 2 F 1 (a,b,c;x) + C (x - 1) 2 Fi(o + 1, b + 1, c + 1; x)] , (28) 
where e is the sign of m and C — ^~j^~r^- Using the transformation 

2 F 1 (a,b,c;x) = rlc-a)r(c-b) 2F t (a,b,a + b - c + 1; 1 - x) 

+ (l-r°' t r( Sy 2 F 1 (c-a,c-b,c-a-b+l,l-x), (29) 
and the x = or x — 1 ~ — 2/r at spatial infinity, we can write the asymptotical form of Y± at spatial infinity as 

Y+ - ie (l> + r^ + D_r~\ m ^j , (30) 
F_ ~i) +r H + £)_ r -IH ) (31) 
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D + = (-2)-IH r ( c ) r ( c - a -_ 6 ), D =( - 2 )H r ( c ) r ( a + ^- c ), (32 ) 

r(c-a)r(c-b) r(5)r(h) 



where 



This matches the asymptotical behavior already obtained in Eq. 

When \m\ < Va, it is easy to see that the boundary condition r~ 1 Y± — » at spatial infinity can be satisfied 
automatically, thus we have no limit from the boundary condition requirement on the value of u>. This is quite 
unusual since the QNM with any frequency, even with the positive imaginary part can be allowed. This implies that 
when \m\ < the perturbation may have growing modes and the black string is unstable against the Fermonic 
perturbation. _ 

On the other hand for the perturbations of \m\ > -y/A, the fulfillment of boundary condition r~ 1 Y± — > at spatial 
infinity calls for the condition that D + = 0. Thus we need to let c — a be zero or negative integer, which gives the 
discrete values of QNM frequencies 

uj = -i(l/2 + n+ \m\/Vxj , n= 0,1,2,... (33) 

where we have recovered A. The purely negative imaginary frequencies indicate that when |m| > VA the perturbation 
of massless Dirac field exhibits the decay mode showing that the black string is stable under such perturbation. 

III. GRAVITATIONAL PERTURBATION 

Now we consider the linear perturbation of the BTZ black string spacetime, which we will denote as gAB + h-AB 
where gAB stands for the components of the unperturbed black string metric and Hab is the metric perturbation. 

A. Perturbation equations 

The perturbation equations can be greatly simplified by choosing appropriate gau ge c onditions, which was discussed 
in detail in studying the gravitational wave in asymptotically flat spacetimes (e.g., [lg, §4.4). We will employ gauge 
conditions on Hab including the tracelessness condition, the transversality condition, namely h\ = 0, h AE = 0, and 
hit* = 0, so that 



h AB (t,r,ip., z) = e nt ~^ 



/ H tt (r,z) H tr (r,z) H tlp (r,z) 

H rt (r,z) H rr (r,z) H rip (r,z) 

H vt (r,z) H^rlr.z) H vv (r,z) 

V H„(r,z) 



(34) 



The linearized Ricci tensor has the form SRbd = —^O^'hsD + RE{Bh^ — RBCDEh CE . Applying this to linearize 
Eq.©, we get 

D^h BD + 2R BCD Eh CE = 0. (35) 

Inserting the perturbation Eq. (|34l) into the above equation, we find immediately that with the trace- free condition, the 
equation governing H zz is decoupled from the rest, giving the description of scalar part of the perturbation. Taking 
the ansatz H zz (r, z) = H zz (r)u m (z), we have 

(Xr 2 - k)H zz . rr + r-^SXr 2 - k)H zz . r - [(Ar 2 - fc)- 1 ^ 2 - 2A + /i 2 r~ 2 + rn 2 ]H zz = (36) 

where m 2 is the effective mass of the Kaluza-Klein modes satisfying 

a^ 1 {a! i u m ) 1 — 6(a) 2 u m + m 2 u m = 0, (37) 

where the prime denotes the derivative with respect to z. In the vicinity of the horizon and near the spatial infinity, 
solutions of Eq. (I36p asymptotically have behaviors near the horizon and the spatial infinity 

H zz ~ (r-r H f^, r^r H , (38) 



-i±y/ m yx-i^ ^oo. (39) 
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Vh = yJk/X denotes the event horizon. In order to have only ingoing wave near the horizon, we require the sign in 
the exponential of Eq. (|38| to be positive. 

Now we turn to discuss the tensor perturbation. Since H zz is independent of the rest, i.e., of the tensor perturbation, 
we set H zz = 0. However to obtain the perturbation equation, we need to take further simplification. Considering 
just the s-wave [l[ that H tip = and H rtp = 0, we have 



H tt (r,z) H tr (r,z) 
h AB (t,r,z) = e™\ H ^ Hrr(r,z) 




(40) 



Using the metric Eqs.Q© an d the perturbation in the form (|40|) . we find that the perturbation equation reduces 



to 



a- 2 (p^h^ + 2Rf Xvp h x ^ + a (a,-% u )' - la^d'h^ = 0, 



(41) 



where the prime denotes the derivative with respect to z, and the greek index relates to the three-dimensional 
spacetime Eq.([5]). Separating the perturbation into /^(t, r, z) — h^,y(t 7 r)v m (z) — e nt H PjV (r)v m (z), we can rewrite 
(l4"T) into 



a (a 1 v' m ) — 2a 1 a"v m + a 



2 2 

m v r , 



0, 



and 



n (8) V + 2R% p h xp - m 2 h^ = 0. 



With the aid of gauge conditions, we can obtain a second order differential equation of H tr from Eq. 
tedious calculations, which reads 



(42) 

(43) 
through 



[n 2 + (Xr 2 - k)m 2 + X 2 r 2 - 2Xk] H t , 



7Xr 2 - k 
r(Xr 2 — k) 



n 2 



5Ar 2 - k 



A(5A 2 r 4 - 13Afcr 2 + 2fc 2 ) 



r(Ar 2 — k) 



Htr,r 



A(3Ar 2 - 2k)(X 2 r 4 - AXkr 2 - k 2 ) 2X 2 r 4 - 2Xkr 2 - k 2 



r 2 (Ar 2 — k) 2 



rO 2 - 



n 2 



Xr 2 



r 2 (Xr 2 — k) 



H t r = 0. 



(44) 



r 2 (Xr 2 -k) 2 ' 

Its asymptotic behaviors near the horizon and the spatial infinity have forms as 

H tr ~(r~r H ) r — *r H , 
H tr ~r- 2± ^ m2 ' x +\ r^oo. 

Considering only the ingoing wave near the horizon, we will just keep the positive sign in the exponent of Eq. 

We have obtained equations describing the scalar part of the gravitational perturbation (]36[) and the tensor pertur- 
bation (|44[) . In the following we will adopt the boundary condition that at the spatial infinity the wave functions H zz 
and H tr vanish so as to solve these equations and examine the stability of the BTZ black string against gravitational 
perturbations. 



(45) 
(46) 



B. Exact Solution of Scalar Part Perturbation 



Taking H zz (r) = r 1 l 2 i\)(r) and to = ifl, Eq. (f3"6"l) can be put in the form of a Schrodinger equation, 

where r* is the tortoise coordinate defined by = J drj (Xr 2 — k)= 2 ^ xk In ^ , such that the horizon is at 
and the spacial infinity locates at =0. The effective potential reads 



(47) 



- DC 



V(r) = (Xr 2 - k)[m 2 - 5A/4 + (fc/4 + ^ 2 )r~ 



(48) 
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For simplicity we will take A = 1. The potential is very similar to that of the massless scalar field perturbation described 
by the Klein-Gorden equation in the BTZ black hole on the brane where V(r) = (\r 2 — k){m 2 — 5A/4+(fc/4+u 2 )r~ 2 ] [20l]. 
which vanishes at the horizon but diverges at spatial infinity. 

Adopting x = l/cosh 2 r% = (r 2 - r^)/r 2 [H HJ, and tp(r(x)) = (1 - xfl^^^-^x- 1 ^ l 2 y{x), we can transform 
Eq. (|47|) into a hypergeometric equation as Eq. (|26|) , but now with a = | (1 — \Jm 2 — \ — iui + ifi) , b = | (1 — \Jm 2 — f — 
io; — i/i), and c = 1 — iuj. To meet the boundary condition at the horizon x — where there are only ingoing waves, 
we choose the solution of the hypergeometric equation to be simply y(x) =2 F\{a, 6, c; x), so that 

i/j = (1 - s )V4-VSP=T/a a .-i«/a 2Fi(fl; 6j q ^ (49) 

il„ = (f - I )V2-Vm a =T/a a .-ia. /2 aFl(0) ft) c . (50) 

We will check the restrictions of the boundary condition at spatial infinity. We write the asymptotical form of H zz at 
spatial infinity in the standard form Eq. (|39[) . Using the transformation relation of hypergeometric function Eq. (I29[) 
we have, at spatial infinity, 



H zz ~ G+r- 1 ^" 1 "- 1 + G_r- Wm2 -\ (51) 



where 



_ T(c)T(c-a-b) _ T(c)r(a + b-c) 

U+ ~ T{c-a)T(c-b)' r(a)r(6) ' { ' 

When m 2 < 2A, it is easy to see that the boundary condition that H zz vanish at infinity is satisfied automatically, 
which is similar to what we observed in the fermonic perturbation. This means that there is no limit on the frequency 
of the perturbation, which can allow the perturbation with growing modes to exist and make the spacetime to be 
unstable. 

When m 2 > 2A, we need G + = to ensure the vanishing of H zz at spatial infinity. Thus we should have c — a or 
c — b to be zero or negative integers, which leads to 



w = ±H- i(2n + 1 + y/m 2 /\ - 1), (53) 

where we have restored the A. Notice that the imaginary part of the frequency is negative, which describes the decay 
mode of the perturbation. Thus the BTZ black string is stable in this case. 

C. Numerical Result of the Tensor Perturbation 



Now we try to solve the equation of the tensor perturbation Eq (|44|) . It has to satisfy the boundary condition that 
Htr vanishes at spatial infinity. 

From the asymptotic behavior Eq. (|46[) . it is easy to see that the boundary condition can be satisfied automatically 
when m 2 < 3A. This again means that in this situation, there is no restriction on the frequency of the perturbation 
and even the growing mode is permitted which can make the spacetime unstable. 

When m 2 > 3 A, we need to solve Eq. ([4"4")l to examine the perturbation behavior. This equation cannot be solved 
analytically, which requires numerical calculation. Considering that the spacetime has the AdS property, we will 
employ the numerical technique developed in [I?) and applied in many AdS black hole spacetimes [lS|. The radial 
equation Eq (|44[) keeps the same form under the transformation r — > ar, m — > raj a and Q — > Cl/a. We will take A = 1 
for simplicity. Taking x — r" 1 , we can rewrite Eq. (|44| in the form 



t(x) u(x) 

Htr.x + 7 \ 

X H [X-XH) 



s(x)H tr , xx + -±J-H tr}W + - \' H tr = 0, (54) 



where xh = ^/rn = 1, and 



s{x) = (2 + rn 2 -0 2 )a; 6 + (4 + 2m 2 -2^ 2 )a; 5 -2(l + m 2 )x 3 -(H-m 2 )a; 2 , (55) 

t{x) = (2 + rn 2 -^ 2 )(a; 6 +a: 5 )-(7 + 2TO 2 -50 2 )(a; 4 +x 3 )-3(l + m 2 )(a; 2 +x), (56) 
u(x) = -(2 + m 2 -fl 2 )x 6 -(5 + m 2 -m 4 + 2m¥ - VL A )x A 

+2(7 + 2rn 2 - m 4 - 3fl 2 + m 2 n 2 )x 2 - 3 - 2m 2 + m 4 . (57) 
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We can expand these functions around xh, e.g. s(x) = Y^=o s n( x ~ x if) n , and so does t(x) and u(x). H tr can also 
be written in the expansion form H tr — Y^=o a n( x ~ x H) n+p - Inserting all these expansions into Eq. (j5~4")) . we can 
solve the leading order equation and find p = — 1 ±fi/2. We will take the plus sign which stands for the in-going wave 
near the horizon. Equating the coefficients of {x — xu) n , we have the recursion relation for a n : 

^ n— 1 

a n = - — ^2[(k + p){k + p - l)s„_jfe + (k + p)tn-k + u n - k ]a k , (58) 
Fn k=o 

where 

P n = (n + p)(n + p - l)s + (n + p)t + u . (59) 

Considering the boundary condition that H tr \ x =o — at infinity, we require that Y^Lo a n(^> m ){~ x H) n — 0, which 
is satisfied only for discrete values of O. In the later discussion, we will take u> = We need to find the zeros of 
J2^Lo a n( tJ 7 TO )( — x h) h = in the complex uj plane. This is done by truncating the series after a large number of 
terms and computing the partial sum as a function of u). We can then find zeros of this partial sum, and check the 
accuracy by seeing how much the location of the zero changes as we go to higher partial sums. 

To check the numerical scheme, we will first use it to calculate the scalar part perturbation, which can be compared 
with analytical results. For the scalar part, the corresponding s(x) = 4x 4 + 8a; 3 + 4a; 2 , t(x) = 8a; 4 + 8a; 3 , and 
u(x) = (1 + 4^i 2 )a; 4 - (6 - 4m 2 + 4/i 2 + 4Sl 2 )a; 2 + 5 - 4m 2 . We truncate the infinite series at N = 200 and get 
the numerical result listed in Table 1. We find that the numerical calculation agrees well with the analytic solution, 
especially for the lower modes. This gives us confidence to use the numerical scheme to explore the tensor perturbation. 
Results for the tensor perturbation are shown in Table 2. It is interesting to find that for m 2 > 3A, all modes have 
negative imaginary parts of QNM frequencies, which shows that in this case the spacetime is stable against tensor 
perturbation. Furthermore we observed that for lower modes the QNM frequency only has the imaginary part in the 
tensor perturbation and surprisingly its value satisfies uj — —i(2n — I + \J m 2 + 1), n = 0, 1, 2.... 





M 


n 


^num 


^cxact 


2 








-2.0000i 


-2i 






I 


-4.0001« 


-4* 






2 


-6.0002i 


-6* 






3 


-8.0003i 


-8* 


2 


1 





±1.0000 - 2.0000« 


±1-2* 






I 


±1.0000 - 4.0000« 


±1-4* 






2 


±0.97- 5.997z 


±1-6* 


3 








-2.4140i 


-2.41421*' 






1 


-4.4142* 


-4.41421* 






2 


-6.41* 


-6.41421* 



Table 1. Comparison between the numerical result and the analytic result of u) for the scalar part perturbation. 



2 

m 










-1 + Vm 2 + 1 


3.001 


-1.0002500* 


-3.0002500i 


±2 x lO" 7 - 


5.0002500i 


1.0002500 


4 


-1.236068* 


-3.236068i 


±0.00007 


- 5.236068i 


1.2360680 


5.25 


-1.5000000* 


-3.5000000i 


±0.00004 - 


5.5000000i 


1.5 


6.25 


-1.6925824* 


-3.6925824i 


±0.00001 - 


5.6925824i 


1.6925824 


8.01 


-2.0016662i 


-4.0016662i 


±3 x 10~ 8 - 


6.0016662i 


2.0016662 


11.25 


-2.5000000i 


-4.5000000i 


±0.000001 - 


6.5000000i 


2.5 


16 


-3.1231056i 


-5.1231056*' 


±4 x 10~ 8 - 


7.1231056*' 


3.1231056 


19.25 


-3.5000000i 


-5.5000000i 


±4 x 10~ 8 - 


7.5000000i 


3.5 



Table 2. Quasinormal frequencies of tensor perturbation obtained through numerical calculation when m 2 > 3A. 
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IV. CONCLUSIONS AND DISCUSSIONS 

We have studied the stability of the BTZ black string against fermonic perturbation and gravitational perturbation. 
It is interesting to see that the BTZ black string can be dynamically stable provided that the m 2 , which is determined 
by the compactification of the extra dimension, is over a threshold value, namely A, 2A, 3A respectively for fermonic 
perturbation, scalar part of the gravitational perturbation and the tensor perturbation. When m 2 is smaller than 
this threshold value, the perturbation can have a growing mode, which indicates that the BTZ black string can be 
unstable and pinch-off to form a black hole. 

The situation that the stability and instability coexist in the dynamical system for the BTZ black string can be 
understood from thermodynamical arguments, if we compare the black string phase with the localized black hole 
phase with fhA — ► 0. The entropy argument provided in [10( shows that the BTZ black string is not a privileged 
stable phase. Its entropy equals to that of the localized BTZ black hole with fhA — » 0. Our dynamical study supports 
this result. 

It is of great interest to extend our discussion to the solution with localized black holes on the brane [10( where the 
geometry induced on the two-brane is described by ds 2 = — (Xr 2 — k — 2fhA/r)dt 2 + (Ar 2 — k — 2fhA/r)~ 1 dr 2 + r 2 d(p 2 . 
The extra term 2fhA/r came from the four-dimensional nature of the black hole. From the entropy argument it was 
argued that when fhA is deviated from zero, entropies of black string and black hole will be nonzero and the black 
string entropy will exceed that of the black hole for small fhA. But for large fhA, the black string entropy will be 
lower. Thus the black string will first be stable and then become unstable and pinch-off to form a black hole [Icj |. 
The dynamical description of this process needs careful investigation and we will report its result elsewhere. 
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